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OPTIMAL CRUISE TRAJECTORIES FOR SUPERSONIC AIRPLANES 
by Fred Teren and Carl J. Daniele 
Lewis Research Center 

SUMMARY 

Equations are derived for maximizing range for specified initial and final values of 
mass and altitude. Constant -velocity flight is assumed, and normal acceleration is 
neglected. 

The optimization problem is solved by using the maximum principle. The resulting 
optimal control histories are shown to consist at most of one maximum, one minimum, 
and one interior thrust arc. Optimal trajectories are obtained, and uniqueness is dem- 
onstrated. Results are obtained for a supersonic airplane. Curves are presented which 
can be used to obtain the optimal trajectory and maximum range for a range of initial and 
final mass and altitude. In addition, the optimal range is compared to the range obtained 
by using the standard cruise trajectory profile (consisting of a Breguet cruise plus 
maximum- and minimum -thrust connecting segments) and to the range at constant- 
altitude cruise. The difference between optimal and standard trajectory range is shown 
to be small. 


INTRODUCTION 

The optimization of airplane trajectories is generally considered as three separate 
optimization problems - takeoff and climb, cruise, and descent and landing. The optimal 
climb problem has been treated in the literature by using the energy -state approximation 
(ref. 1), as well as by treating a more precise formulation of the system equations with 
the attendant complexities (ref. 2). The descent phase is generally highly constrained, 
and therefore not subject to mathematical optimization. The importance of optimizing 
the cruise trajectory (maximizing cruise range) is obvious. Yet, this phase has re- 
ceived little attention in the literature. The reason for the lack of attention is that the 
approximations generally made in considering the cruise trajectory optimization reduce 
the solution to a nearly universal equation for all airplanes, the so-called Breguet range 



equation. Furthermore, the Breguet cruise seems intuitively to be difficult to improve 
on, since the associated approximations seem to be quite reasonable. 

The Breguet range equation is derived by assuming steady flight. All accelerations 
are assumed to be zero, and the climb rate is also neglected. Therefore, altitude be- 
comes a control variable for the problem, and the Breguet range equation results. In 
the present study, the cruise trajectory model is made more realistic by including the 
altitude as a state variable. With this formulation, the cruise trajectory is defined as 
the portion of the flight at cruise velocity, with specified initial and final altitudes corre- 
sponding to the end of climb and the beginning of descent. The optimization problem is 
formulated by using the maximum principle (ref. 3). Costate equations are derived and 
are used to determine the optimum thrust level as a function of the state variables. 

The procedure for obtaining optimal trajectories is illustrated for a supersonic air- 
plane. A set of curves is presented from which maximum range can be determined for 
any values of initial and final mass and altitude. Comparisons are also made between 
the optimal range and the range obtained by using two suboptimal cruise trajectory 
profiles . 


CRUISE RANGE OPTIMIZATION 
Equations of Motion 

It is desired to maximize range for an airplane with fixed initial and final mass. 
The velocity is assumed to be constant, and the initial and final altitudes are specified. 
The flight time is free. The equations of motion are 


v = — (T - D - mg sin y) 

(la) 

m 

y = — — (L - mg cos y) 

(lb) 

mv 

h = v sin y 

(lc) 

m = - — S 

(Id) 


g 

All symbols are defined in the appendix; some of the variables are illustrated in figure 1. 
The thrust is assumed to be parallel to the airplane longitudinal axis , and the angle of at- 
tack is neglected. The specific fuel consumption S is a function of thrust and altitude. 
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Equations (1) apply for planar flight with a constant gravity field (that is, a flat Earth is 
assumed). The range is given by 


R = / * v cos y dt (2) 

*0 

Breguet Cruise 

The Breguet equation has generally been used to calculate airplane cruise range. 

For completeness and for comparative purposes, the Breguet equation is derived in this 
section. The following assumptions are made: 

(1) The normal acceleration vy is negligible. This does not mean that y is con- 
stant, only that the acceleration required to change y is negligibly small. 

(2) The vertical velocity h is negligible. Here again, h is not necessarily con- 
stant. Neglecting h requires setting sin y = 0 and cos y = 1 in equation (1). Further- 
more, this assumption makes h a control variable. 

Using these assumptions, along with constant velocity, gives for equations (la) 
and (lb) 


T = D 

(3a) 

L = mg 

(3b) 


and the range becomes 


R = 


A f r m f 

j vdt= -y 

* / m n 


vg dm 
TS 


c‘Y = - 


TS 


This equation can be rewritten by using equation (3) as 


R = - 



( 4 ) 


which is the Breguet range equation. Maximum range for the Breguet cruise is obtained 
by maximizing, at each flight time (or mass), the value of L/DS. The maximization 
may be performed independently at each point along the trajectory since altitude is as- 
sumed to be a control variable and may be changed discontinuously if required. 
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At this point, it is useful to describe the airplane and atmospheric characteristics in 
more detail. The following assumptions are made: 

(1) A parabolic drag polar is assumed; that is, 

c d = c d,o + kc l 

where C D Q and K are constants, independent of altitude. The lift and drag are given 

by 


L = C l QA 
D = C d QA 

where A is the reference area and Q is the dynamic pressure, defined by 



(2) An isothermal atmosphere is assumed; that is, 


(5a) 

(5b) 


(6) 


P = PQ e 


-/3h 


(V) 


This also implies that speed of sound is constant. Hence, Mach number remains con- 
stant during a constant -velocity cruise. 

(3) The airplane specific fuel consumption is a function (for constant -velocity flight) 
only of the thrust coefficient Cp, where 


C 


F 


T 

QA 


(8) 


In addition, Cp has upper and lower limits, Cp u and Cp ^ , which are independent of 
altitude. 

With these assumptions, the operating conditions for maximizing Breguet cruise 
range can be determined more explicitly: 


L _ C l^ A _ C L 
DS C d QAS C^S 


4 


Equations (3a), (5b), and (8) can be combined to give Cp = Cp. Also, Cp can be ex- 
pressed in terms of Cp as follows: 


C 


L ~ 




( 9 ) 


Therefore, 


L _ 1 ^/ C F C D,0 

DS C F S » K 


and the optimum Cp can be obtained by differentiating L/DS with respect to Cp 
(with S = S(Cp)) and setting it equal to zero. This results in the following equation, 
which can be solved for the optimum value of Cp: 


c F s 

2 


(Cp 


- c 


D,0 J 


+ c 


F 



(10) 


It should be noted that the optimum value of Cp, is independent of mass, and hence con- 
stant for the entire cruise trajectory. Also, from equation (9), Cp is constant; thus, it 
follows that L/D is also constant. In addition, since Cp = L/QA = mg/QA, 
m/Q = ACp/g = constant. Thus, Q must decrease (altitude increases) as mass de- 
creases. The Breguet cruise range can be obtained by integration of equation (4) and 
using the fact that L/D is constant: 


R = v 




(ID 


Inclusion of Altitude Effect 

In the derivation of the Breguet range equation, altitude was assumed to be a control 
variable and was allowed to change discontinuously if required. At this point, the prob- 
lem will be stated more realistically by including altitude as a state variable. The state 
equations (1) may be rewritten by using equations (5) to (8) as 
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(12a) 




EE sin y 
AQ ) 


y = — (C l QA - mg cos y) 
mv 


(12b) 


2Q 2 A 


Q = -/3Qv sin y + (c F - C D Q - KC 2 


mv 


El sin y\ 
QA / 


(12c) 


mz.^CjSfCp) (12d) 

g 

where the state variable h has been replaced by the equivalent state variable Q. 

In order to simplify the problem, it will be assumed here (as it was in the derivation 
of the Breguet range equation) that normal acceleration is negligible, cruise velocity is 
constant, and y is small (i. e. , cos y « 1). However, the state variables Q and m 
will be retained. Setting equations (12a) and (12b) equal to zero, solving for sin y, and 
\asing this result in equation (12c) yield 


__ (3Q^vA L 


Q= - 


Cp - C 


Kg 2 m 2 


mg 


\ F D >° 


Q" 


The optimization problem can now be stated as 


maximize /. v dt 
r 0 


f- 

'Jt n 


subject to the differential constraints 


Q = - ffiiWc. - c D 0 - 

"! \ F D ’° aV 


m = -^C F S 


(13) 


(14) 


(15a) 

(15b) 


an inequality constraint on the control Cp 

C F, l ~ C F ~ C F,u 
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and boundary conditions on Q and m. The problem is formulated by using the 
Pontryagin maximum principle in minimum form (ref. 3). The variational Hamiltonian 
can be written as 


H = -v - X. 


/3Q 2 vA 

mg 


' C F " C D, 0 


Kg 2 m 2 ^ 

A 2 Q 2 


- X 


2 — C F S + /i(C F - C FjU )(C F - C F ^) 


(16) 


The costate equations are 


Xj = 2 Xj 


/3QvA 

mg 


(C F - C jj q) + — ac f s 

* O’ 


(17a) 


x 2 - - 


X^QvA Kg 2 m 2 ^ 

— I c f-°d,o + ~TT 

gm^ \ A Z Q Z 


(17b) 


The optimal C F is obtained from 


8H 

dC-rp 


= -X £9?^ -X 2 AQ-1- (SC F ) + m(2C f - C F U -C F t ) =0 
1 mg z o- rfp _ r r r f u r 


g dc 


(18) 


and subject to 


8 2 H 

sc| 


— ? AQ (SC-p) + 2/i > 0 

g dC 2 


(19a) 


and 


H = 0 (19b) 

on interior arcs (i. e. , C F ^ < C F < C F u ). °P^ ma l thrust coefficient history may 

consist of maximum, minimum, and interior arcs. It is instructive to determine the 
variation of C F on optimal interior arcs. It will be shown that an optimal interior arc 
is uniquely determined by initial values of C L and C F . On interior arcs , equa- 
tions (18) and (19b) give 
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-X../3QV d 

— X 9 -^_(SC Tr ) = 0 (20a) 

m * dC F * 

For notational simplicity, let a(Cp) = C F S(C F ) and let ' denote differentiation with 
respect to C F . Then equation (20a) may be written 


a’ _ X 1 Q 
j3v ^2 m 

Differentiation of equation (20b) with respect to time gives 


(20b) 


ct c _ u 


d/ X l Q ^ 


Q 


I3v 


dtyXgny Xgm 


2X PQvA (c ^ _ c ^ ft) + J Aa 


mg 


'F “ D, (r 


X 1 /3A r\2. 


X 2 m g 


r. 2 2\ X 1 Q 2 Aa X 2 Q 3 0vA / „ 2 2) 

(!\|c f -c D(| .Sjl.J: 1Z_C F -C n n + Km £ 

2. \ F D ’° . 2„2 / , „2_ ,2„3_ \ F D ’° a 2q2 


A Q / Xgm g X^m g 


Substituting for X^/Xg from equation (20b) results in 


a " CF . - 3. LA (a- .Jc F -C D0+ «2W-\) 

(3v mg /3v \ * U ’ U .2-2 a.’/ 


A 2 Q 2 a 


or, since L = C L QA = mg, 


6 F = ^- ( ^- a X C F- C D,0 + KC L-^ 

L 


(21) 


Thus, the rate of change of C F on an interior arc is uniquely determined by the values 
of C F and C F (for given values of the airplane parameters). The rate of change of 
C L on an interior arc can also be obtained by differentiation: 


C T = -l(^A = -a + /3v( C, 
L dt \QA J 


C D,0- KC L)=K C F- C D,0- KC i-^' 


( 22 ) 


so that the rate of change of C L is also determined by the instantaneous values of C F 
and C L . A constant control arc (C p constant) can exist only if C L is also constant. 
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Thus, 


and 


' C D, 0 + 


KC( 



(23a) 




(23b) 


must be satisfied simultaneously. The value of Cp on this arc is determined by 
solving 


C 


F 


" C D, 0 



The value of C ^ can then be determined from equation (23b) as 



(24a) 


(24b) 


It is interesting to compare the constant control arc given by equation (24) with the 
Breguet cruise equations (9) and (10). If equation (10) is rewritten in terms of a, the 
result is 


C 


F 


~ C D,0 



(25) 


The similarity of equations (24) and (25) is evident. In fact, equation (24a) represents a 
modified Breguet cruise - modified by the additional thrust required to change altitude. 


Construction of Cp - Cp Diagram 

Because of the nonlinearity of the variational equations required for maximizing 
range, the optimal trajectories are very difficult to obtain by conventional methods 
(e. g. , solution of the two-point boundary value problem). It was found that the important 
features of the optimal trajectories could be conveniently displayed on a graph of lift co- 
efficient as a function of thrust coefficient. From these features, an approximate (but 


9 



quite accurate) procedure was derived for obtaining the optimal trajectories. In the fol- 
lowing development, the Cp - diagram is presented and discussed, and the solution 
procedure is derived. 

Equations (21) and (22) can be used, for a given airplane-engine configuration, to 
construct a Cp - Cp diagram of all possible optimal interior arcs. The airplane- 
engine data used here are representative of a supersonic transport (SST). The assumed 
values of the airplane parameters are given in table I, along with the inverse atmo- 
spheric scale height j3. 

The specific fuel consumption S is shown as a function of Cp in figure 2. The 
value of S is constant until afterburning is initiated at Cp, after which S increases 
with Cp. The data in table I and figure 2 apply at a velocity of 754 . 3 meters per second 
(corresponding to a Mach number of 2.7). 

From equations (21) and (22), define 

f(Cp) = Cp - Cj-j q — - (26a) 

b(Cp) = Cp - Cpv q (26b) 

r r J ’ V Rv 


Then equations (21) and (22) become 

C F — (/3v - a’)(f + KCx 
a"C L V L 

C L = pv(b - Kcty 


(27a) 


(27b) 


Also, 


t T 

f T aa 
(a’) 2 


(28a) 


b’ = 1 - — (28b) 

/3v 

Since a" > 0 and a' < /3v for the values assumed, f and b are monotonically in- 
creasing functions of Cp (see eqs. (28)). Also, f(Cp) is negative for all values of Cp, 
while b(Cp ^) is negative and b(Cp u ) is positive. These characteristics define the 
shape of the Cp = 0 and C L = 0 curves, as shown in figure 3. Above the Cp = 0 
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curve, C F increases; below, Cp decreases. Above the C L =0 curve, Cp de- 
creases; below, Cp increases. The intersection of the Cp = 0 and Cp = 0 curves is 

the modified -Breguet -cruise (MBC) point, as discussed previously in connection with 

• • 

equations (24) and (25). For each Cp, Cp combination, the values of Cp and Cp are 
uniquely determined (although not shown explicitly). Thus, for a given Cp, Cp start- 
ing point, the entire optimal trajectory is determined from that point on. 

In figure 3, a discontinuity in Cp exists between Cp and Cp This discontinu- 
ity is due to the constant S assumed for values of Cp less than Cp, which results in 
a" = 0 for these values and makes Cp increase to infinity in equation (27a). Thus, 
when Cp decreases to Cp on an interior arc, it will jump to Cp Alternatively, 
trajectories starting at Cp ^ will jump to Cp if an interior arc is required. 

Four regions are defined in figure 3 by the Cp = 0 and Cp = 0 curves and their 
intersection. In region I, Cp decreases and Cp increases; in region n, Cp and Cp 
both decrease; in region III, Cp increases and Cp decreases; and in region IV, Cp 
and Cp both increase. The following optimal trajectories are illustrated in the figure: 

(1) Trajectory A starts in region I. First Cp decreases while Cp increases to 
the Cp = 0 curve. Then Cp decreases and Cp decreases to Cp; Cp then jumps to 
Cp j and Cp decreases. 

(2) Traj ectory B starts in region I. First Cp increases while Cp decreases to 
Cp = 0. Then Cp increases while Cp increases to Cp . Finally, Cp increases 
while Cp remains at Cp . 

(3) Trajectory C starts in region HE. As Cp decreases, Cp increases to Cp = 0. 

Then Cp decreases while Cp decreases to Cp ^ and remains at that value. 

(4) Traj ectory D starts in region HI. As Cp increases, Cp decreases to the 

Cp = 0 curve. Then Cp increases and Cp increases to Cp . Finally, Cp in- 

creases while Cp remains at Cp u . 

For all trajectories starting in region II, Cp decreases and Cp decreases to Cp, 
while in region IV, Cp increases and Cp increases to Cp . 

An additional question that remains is the number of interior arcs which can occur 
for a single optimal trajectory. A related question is whether the value of Cp can jump 
discontinuously from Cp ^ to Cp u or vice versa, or if interior arcs must be used 
for this transition. 

The former question can be answered by referring to figure 3. From this figure it 
can be seen that all interior arcs which terminate at C--, „ arrive at C™ in region IV. 

„ . r jU r 9 u 

In this region, Cp is positive, so that the trajectory can never return to region I (along 
the line Cp = Cp ) to start another interior arc. Similar reasoning shows that once 
Cp has reached Cp ^ (regions I or II), it can never return to region in to start another 
interior arc. 
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Additional interior arcs might still exist if Cp could jump dis continuously from 
Cp j to Cp u or vice versa. Then, for example, Cp could jump from Cp ^ to 
Cp while Cp = 0. 01, and a second interior arc could start at that point. To show that 
this is not possible, consider first equation (16). The control Cp is chosen to minimize 
H, and if Xg is positive, H can always be made negative by choosing Cp such that the 
coefficient of Xj has the same sign as Xp However, H = 0 for a free-time optimum 
solution. Therefore, Xg must be negative. Now if Cp can jump from Cp u to Cp j, 
H(C F u ) = H(Cp j) at the transition time. But from equation (19), 



X 9 AQa 

— + 2ju > 0 

g 


so that H(Cp) is convex. Therefore, if H(Cp u ) = H(Cp ^ ), H(Cp) < H(Cp u ) for every 
intermediate value of Cp, as shown in sketch (a). 



(a) 


Hence, Cp would be equal to Cp q, and not Cp u or Cp j. Therefore, no jumps 
from Cp u to Cp ^ (or vice versa) are possible. Furthermore, the convexity of 
H(Cp) verifies the optimality of interior arcs, since it demonstrates that the necessary 
condition in equation (19) is satisfied for such arcs. 

The previous discussion proves that the control history for an optimal trajectory 
consists (at most) of one maximum, one minimum, and one interior arc. Furthermore, 
the optimal trajectory is unique. 
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RESULTS AND DISCUSSION 
Long-Duration Cruise 

As noted earlier, the values of Cp and do not change with time at the MBC 
point. Although not shown explicitly in figure 3, the values of Cp and Cp change 
slowly with time near the MBC point, but change quite rapidly away from this point. For 
example, trajectory A in figure 3 requires 6. 25 minutes and has a range of 283 kilo- 
meters. The initial and final altitudes are both 9. 15 kilometers. Trajectory B requires 
only 1. 54 minutes to climb from 9. 15 kilometers to Cp at an 18. 29 -kilometer alti- 
tude. The range traversed is 69.4 kilometers. Trajectory C starts at 21. 37 kilometers 
and requires 8.42 minutes to descend to 9. 15 kilometers, while covering a range of 
381.4 kilometers. 

While these trajectories are optimal, they are not normally of interest since an SST 
must cruise for long periods of time (several hours). Since trajectories A to C have 
durations of several minutes, it is clear that optimal cruise trajectories of several hours 
duration must spend much time very near the MBC point. In fact, it is reasonable to as- 
sume that such trajectories include a long -duration phase at the MBC point, together 
with connecting arcs to the initial and final conditions. 

Figure 3 may be used to determine how such trajectories can be numerically con- 
structed. A comparison of trajectories A and B in region I indicates that there is a 
unique trajectory, between A and B, which terminates at the MBC point. All other tra- 
jectories in region I diverge from this point. Similarly, a single trajectory exists in 
region III which terminates at the MBC point, which lies between trajectories C and D. 
Also, there are only two trajectories which leave the MBC point, and these lie in regions 
n and IV. The four trajectories of this type are shown in figure 4. The optimal trajec- 
tories then start with connecting arc 1 or 3 (depending on the initial altitude), continue 
with an MBC, and conclude with connecting arc 2 or 4, again depending on the desired 
final altitude. Since trajectories generally consist of climb, cruise, and descent, con- 
necting arcs 1 and 2 will normally be used. Connecting arc 1 is the climb -transition arc 
between the starting point and the MBC point, while connecting arc 2 is the descent- 
transition arc. These transition arcs correspond to those phases of the trajectory in 
which the airplane is at cruise velocity but not at optimum altitude for the MBC. 

It has been established that the long-duration cruise trajectories of interest consist 
essentially of two definite connecting arcs (with end points depending on the initial and 
final values of mass and altitude) and an MBC. It remains to determine the values of 
Cp, Cp, and m, as functions of time, on connecting arcs 1 and 2. For connecting 
arc 1, these values are best determined by numerical integration backwards in time from 
the MBC point. For connecting arc 2, the integration is performed forward from the 
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MBC point. In both cases, the initial values of C L and Cp are perturbed a small 
amount from the MBC values into the proper region. This procedure gives small, but 
nonzero, values of Cp and Cp which are needed to start the integration of equa- 
tions (21) and (22). 

The resulting values are presented in figures 5 to 10. The initial value of Cp, 

C L q, is shown as a function of time spent on connecting arc 1 in figure 5. The required 
value of Cp q is a function of the initial altitude and mass. For convenience, Cp/m is 
plotted as a function of initial altitude in figure 6. From this figure, Cp q can be easily 
determined. The ratio of the mass at the start of the MBC to the initial mass m g]3(V m 0 
is also a function of connecting arc time t gB £. This may be seen by rewriting equa- 
tion (15b) as 


m = Z^CjS = - 


ma 



Since Cp and Cp are definite functions of time, this equation may be integrated to 
give 


m. 


m, 


SBC 



a(C w ) 


dt 


The result is shown in figure 7. Thus, connecting arc 1 is completely determined for 
given values of initial mass and altitude from figures 5 and 7. Similar data are pre- 
sented for connecting arc 2 in figures 8 and 9. Figure 8 shows Cp f as a function of 
time spent on connecting arc 2, while figure 9 gives the mass ratio m EB( Vm f . 

For given values of initial and final mass and altitude, the time spent on connecting 
arcs 1 and 2 is determined, as are the initial and final mass at the MBC point. It re- 
mains to calculate the time spent on the MBC arc t B( ~,. This time may be calculated by 
integrating equation (15b), which results in 


^C - 



log(^ 

\ m EBC/ 


(29) 


where the values of Cp, Cp, and S are those of the MBC arc. The results of equa- 
tion (29) are plotted in figure 10 for convenience. However, for greater accuracy, equa- 
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tion (29) should be used in calculating tg C . For the parameters assumed in this report, 
(CjVCjS)^ = 4. 8972 hours. 

In order to illustrate the use of the figures, an example is given. Suppose 

m Q = 291 000 kg 
irij = 184 000 kg 
h 0 = hf = 9. 2 km 

Then from figure 6, 

Ct rj < 

— = 1.0X10"' kg" 1 
m 

C L o = 1. 0X10 " 7 m Q = 0.0291 
C L f = 1. 0X10" 7 m f = 0. 0184 

Next, from figure 5, 

tgBc = 10- 35 min 

and from figure 7, 

= 0. 9585 nig — 27 8 924 kg 

From figure 8, 

*EBC = 129 min ^ C L = 

/ S -f 

and from figure 9, 

m EBC = 9 ^02 m f = 193 237 kg 
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Now from equation (28), for m EBc/ m SBC = 692795, 

t BC = 1. 7974 hr 


Therefore, the total trip time is 

t = 10. 35 + 60(1. 7974) + 12.9 = 131. 1 min 
and the total cruise range is 


R = Vt = 5934 km 

It is interesting to compare the maximum range with the corresponding range which 
would be obtained if the standard procedure is used to generate the connecting arcs. With 
the standard procedure, the maximum thrust coefficient is used to climb to the Breguet 
cruise altitude, and minimum C F is used to descend from the Breguet cruise altitude to 
the final altitude. A comparison of standard trajectory range is made with the optimal 
range in figure 11, for m^ = 291 200 kilograms and m^ = 184 200 kilograms. The range 
increase obtained by using the optimal trajectories is plotted as a function of initial alti- 
tude. The initial and final altitudes are assumed to be the same. The results in figure 11 
show that the range improvement is small, less than 7 kilometers even for initial and final 
altitudes as low as 9 kilometers (the optimal initial MBC altitude is about 17. 4 km). 


Short-Duration Cruise 

It is reasonable to expect that the range difference will always be small whenever the 
flight is of long enough duration so that it is advantageous to reach the Breguet cruise 
altitude. However, for very short duration flight starting from a relatively low altitude, 
the entire optimal trajectory is a connecting arc which remains below the Breguet cruise 
altitude, and the range penalty incurred by using the standard trajectory could be 
important. 

Therefore, a comparison of optimal trajectory range with standard trajectory range 
was made for propellant masses ranging from 0 to 17 000 kilograms. The results are 
presented in figure 12, for initial and final altitudes of 9. 15 kilometers and a final mass 
of 184 200 kilograms. For fuel mass up to 6900 kilograms, the standard trajectory con- 
sists only of maximum and minimum thrust arcs, since there is insufficient fuel for the 
standard trajectory to reach the Breguet cruise altitude. The range increment increases 
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to a maximum of about 9 kilometers at a fuel mass of 6900 kilograms. For fuel masses 
greater than this value, the standard trajectory also includes a Breguet cruise segment 
of the required duration. The range increment then decreases rapidly to about 7 kilo- 
meters. The optimum trajectory consists of an interior arc for fuel masses up to 17 000 
kilograms. At this point, the connecting arc passes so close to the MBC point that a 
Breguet cruise phase is included. 

For short-range (low-propellant-mass) cruise, a constant -altitude trajectory might 
reasonably be considered to be a good approximation to the optimum. A comparison of 
optimum cruise range with constant -altitude cruise range is presented in figure 13 as a 
function of fuel mass. As in figure 12, the initial and final altitudes are 9. 15 kilometers, 
and the final mass is 184 200 kilograms. It can be seen from the figure that constant- 
altitude cruise is a poor approximation to the optimum even for small fuel masses. For 
example, the range increment is 510 kilometers for a fuel mass of 15 000 kilograms. 


VERIFICATION OF APPROXIMATIONS 

In the derivation of the optimal trajectories presented, it was assumed that y is 
small and that the normal acceleration is negligible. These approximations may be 
checked at this point by observing the resulting values on the optimal trajectories. 

The value of y may be determined by setting v (eq. (12a)) equal to zero, which re- 
sults in 


sin y = — (C F - C D ) 
mg 

For long -duration cruise, the largest values of y occur on the climb -transition arc. 

The value of y on this arc is shown as a function of time in figure 14 for a trajectory 
starting at an 11. 4 -kilometer altitude. The altitude profile is also presented in this fig- 
ure. For values of y less than 8°, cos y is greater than 0.99, so that the approxima- 
tion cos y = 1 is valid to within 1 percent. Figure 14 shows that y is less than 8° for 
altitudes greater than about 13. 4 kilometers. Furthermore, even if the initial altitude is 
as low as 9 kilometers, the time required to ascend to 13 kilometers is less than 1/2 
minute (out of a total trajectory time of several hours). Therefore, the small angle ap- 
proximation seems to be valid. 

Consider next the approximation of neglecting normal acceleration. In equation (12b), 
y was neglected and was set equal to mg/QA. The effect of the approximation is to 
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neglect the lift required to change y. The true required lift coefficient is given by 


C**C L (cosy+2) 


where 


C 


L 


_ 

QA 


is the approximate lift coefficient. The approximate and true lift coefficients are plotted 
in figure 15 as a function of time on the climb -transition arc. The difference between 
C L and c£ is small near the MBC altitude, but increases to 0. 03 at an altitude of 11. 4 
kilometers. The error in lift coefficient results in an error in drag coefficient (C^ is a 
function of C p), which changes Cp by an equal amount in order to maintain the desired 
climb rate. The true drag coefficient is given by 

c d = c d,o +k < c l> 2 

and is plotted in figure 16 along with the approximate drag coefficient as a function of 
time on the climb -transition arc. Figure 16 shows that the error in C D (hence Cp) is 
fairly small (less than 10 percent for altitudes above 13 km) even though the error in C L 
is large. Furthermore, the amount of time spent in the low -altitude region, where the 
error is substantial, is very short relative to the total cruise time. Therefore, neglec- 
ting normal acceleration represents a reasonable approximation. 


CONCLUDING REMARKS 

The problem of maximizing range for constant -velocity flight between fixed terminal 
conditions (mass and altitude) has been studied. Normal acceleration is neglected, as in 
the derivation of the Breguet cruise, but the thrust required to change altitude is included 
Range is maximized for a supersonic airplane and compared to the range obtained by us- 
ing the standard trajectory, consisting of a Breguet cruise arc and maximum- and 
mi nimum -thrust connecting arcs. The difference between optimum and standard trajec- 
tory range was shown to be small. However, the optimal range can be calculated without 
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great difficulty, as described in this report. Therefore, the method derived in this re- 
port should be considered for use in cruise trajectory calculations. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 13, 1971, 

110-06. 
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APPENDIX -SYMBOLS 


2 

A reference area, m 

a,b functions of Cp introduced for notational simplicity 
Cj-j drag coefficient, nondimensional 

true drag coefficient, nondimensional 
Cd 0 zer °-lift drag coefficient, nondimensional 
Cp thrust coefficient, nondimensional 

Cp minimum Cp for afterburning, nondimensional 

Cp lift coefficient, nondimensional 

Cp true lift coefficient, nondimensional 

D drag, N 

f function of Cp introduced for notational simplicity 

2 

g acceleration due to gravity, m/sec 

H variational Hamiltonian, m/sec 

h altitude, m 

K slope of drag polar, nondimensional 

L lift, N 

m mass, kg 

2 

Q dynamic pressure, N/m 

R range, km 

S specific fuel consumption, sec - * 

T thrust, N 

t time, sec 

v velocity, m/sec 

j3 exponent for isothermal atmosphere, m - * 

y flight path angle, rad or deg 

Xp costate variable, (m )(sec )/kg 

Xg costate variable, m/kg 

p costate variable, m/sec 
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Q 

p density, kg/m 

3 

Pq surfe.ce density, kg/m 

Subscripts: 

BC Breguet cruise 
EBC end Breguet cruise 
f final 

l lower limit 

u upper limit 

SBC start Breguet cruise 
0 initial 

Superscripts: 

derivative with respect to time 
derivative with respect to thrust coefficient 
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TABLE I. - AIRPLANE AND ATMOSPHERIC CONSTANTS 


Velocity, v, km/sec 0.7543 

Reference area, A, m 704 

Slope of drag polar, K 0.5 

Zero -lift drag coefficient, C D Q 0.00878 

Exponent for isothermal atmosphere, km“* 0. 16 



Figure 1. - Definition of trajectory variables. 



Thrust coefficient, Cp 


Figure 2. - Specific fuel consumption as function of 
thrust coefficient. 
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Figure 12. - Range increase for short-duration flight of optimal trajectory over 
standard trajectory as function of fuel mass. Final mass, 184 200 kilogram s f * 
initial and final altitudes, 9.15 kilometers. 



Figure 13. - Range increase for short -duration flight of optimal trajectory over 
constant-altitude (9. 15-km) trajectory as function of fuel mass. 
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Figure 14. - Trajectory variables on climb-transition arc. 
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Figure 15. - Approximate and true lift coefficients on climb-transition arc. 
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